Simultaneous Extrema in the Entropy Production 
for Steady-State Fluid Flow in Parallel Pipes 
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Steady-state flow of an incompressible fluid in parallel pipes can simultaneously satisfy two con- 
tradictory extremum principles in the entropy production, depending on the flow conditions. For 
a constant total flow rate, the flow can satisfy (i) a pipe network minimum entropy production 
(MinEP) principle with respect to the flow rates, and (ii) the maximum entropy production (MaxEP) 
principle of Ziegler and Partridge with respect to the choice of flow regime. The first principle - 
different to but allied to that of Prigogine - arises from the stability of the steady state compared 
to non- steady- state flows; it is proven for isothermal laminar and turbulent flows in parallel pipes 
with a constant power law exponent, but is otherwise invalid. The second principle appears to be 
more fundamental, driving the formation of turbulent flow in single and parallel pipes at higher 
Reynolds numbers. For constant head conditions, the flow can satisfy (i) a modified maximum 
entropy production (MaxEPMod) principle of Zupanovic and co-workers with respect to the flow 
rates, and (ii) an inversion of the Ziegler-Paltridge MaxEP principle with respect to the flow regime. 
The interplay between these principles is demonstrated by examples. 

PACS numbers: 05.70.Ln 47.15.Fe 47.27.Cn 47.27.nf 47.85. Dh 

Keywords: maximum entropy production, minimum entropy production, fluid mechanics, steady state, pipe 
network, dissipative, complex system 



1. INTRODUCTION 

For millennia, the discipline of fluid mechanics has 
inspired many of humanity's deepest thinkers, includ- 
ing Archimedes, Al-Khazini, da Vinci, Torricelli, New- 
ton, Pascal, the Bernoullis, Euler, d'Alembert, Lagrange, 
Laplace, Navier, Stokes, Helmholtz, Reynolds, Taylor, 
Rayleigh, Prandtl, Onsager and Kolmogorov. Despite 
this long and rich history, and the more recent develop- 
ment of thermodynamics over the past 150 years, surpris- 
ingly few workers have sought to establish the connection 
between fluid mechanics and what is generally referred 
to as non- equilibrium thermodynamics: the study of sys- 
tems which experience one or more flows of various phys- 
ical quantities in response to thermodynamic gradients. 
Of particular importance to the latter is the concept of 
entropy production, as defined by [H-Q : 



(i) 



where Fj is the jth thermodynamic force acting on the [§4lQ| 



system and jj is the corresponding flux of the jth phys- 
ical quantity. Examples of the forces include the gradi- 
ents V(l/T), -V(/i c /T), -V(v/T) and E/T, where k 
is Boltzmann's constant, T is absolute temperature, fi c 
is the mass chemical potential of the cth constituent, v 
is the mass-average fluid velocity and E is the electric 
field. These generate, respectively, corresponding fluxes 
of heat j q , mass of each constituent j c , momentum r 



and charge j z 0-0]. The lack of association with fluid 
mechanics is even more surprising with the knowledge 
that non-equilibrium thermodynamics is essentially the 
science of dissipative systems [2|-l5j - those which dissi- 
pate available energy as heat - of which factional fluid 
flow is a prominent example. Remarkably few reference 
books in fluid mechanics even mention the entropy pro- 
duction, whilst many of those that do [e.g. @, 0] examine 
it in passing as merely another transport-related quan- 
tity, seemingly of no great importance. The connections 
between entropy production and the dissipation of avail- 
able energy (exergy), heat generation, work expenditure 
and power consumption are usually not made explicit. 
Part of the problem is that fluid mechanics now tends to 
reside between the disciplines of physics and engineering; 
today, few physics degree programs contain more than 
a single course in fluid mechanics (if any!), whilst few 
engineering programs devote adequate time to "esoteric" 
concepts such as entropy. Of course, exceptions to this 
lack of association do exist. Most relevant to the present 
analysis include several detailed reference works by Bejan 
101, written for a somewhat different purpose (see Sj2]), 
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the quite farsighted work of Paulus Jr and co-authors 
[Til [l2| (which, however, contains some errors), and a 
recent note by Martyushev 

In thermodynamics, undoubtedly the most far- 
reaching advance was the development of extremum prin- 
ciples to determine the equilibrium position of a system. 
These include maximisation of the thermodynamic en- 
tropy of a closed system [l4], [H| and its dual principle, 
minimisation of the free energy of an open system [16U17| . 
In non-equilibrium thermodynamics, several extremum 
methods have been proposed to determine the stationary 
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(steady) state [l8l - [20| . The first was advocated by On- 
sager [2114231 ] , involving maximisation of the entropy pro- 
duction less a dissipative term, widely but inaccurately 
referred to as a "minimum dissipation" method. This 
method is connected to other modified maximum entropy 
production (MaxEPMod) principles for the analysis of 
viscous fluid flows |2-i - t2q ] and electrical circuits |27H30l |. 
The next (and best-known) non-equilibrium extremum 
method involves minimisation of the entropy production 
(MinEP), expounded over many years by Prigogine 
This method has inspired or is connected to other MinEP 
(or minimum heat production) methods for the analysis 
of electrical 0, 13014331 ] and fluid flow systems [ll|, Ell • A 
third, more elusive principle has been proposed in sev- 
eral fields (e.g. by Ziegler and Paltridge) based on em- 
pirical evidence, invol ving maximisation of the entropy 
production (MaxEP) |34l436| . This has been used to 
predict the steady-state behaviour of a variety of non- 
equilibrium "complex systems" , including heat transport 
by convection in heated vessels (Benard cells ) 13711 , plan- 
etary atmospheres [35l - l4fjj and the mantle |4l|: global 
water and nutrient cycles [13, El, EH; crystal growth 
and diffusion [H, EH; resource degradation in ecosys- 
tems [HI, E3 ; biochemical reaction mechanisms [H, EH 
and earthquake frequencies [Eoj]. Other extremum prin- 
ciples have also been given for other purposes, e.g. meth- 
ods based directly on maximising the Shannon entropy 
(MaxEnt) [HHl]- 

The aim of this work is to examine various extremum 
principles for the analysis of fluid flows in pipes, both 
individually and in parallel. The more difficult problem 
of flow in a generalised flow network - in which not all 
nodes need be connected - is deferred to a later analysis. 
Firstly, the various formulations of these principles are 
classified and the principles of fluid flow in a pipe are 
reviewed in light of their connection to frictional losses 
and hence entropy production (Q. This leads into sepa- 
rate analyses of entropy production phenomena in single 
(0 and parallel pipes (©, the latter under both con- 
stant total flow rate or constant head conditions. In a 
single pipe, empirical data provide strong evidence that 
the Ziegler-Paltridge MaxEP principle applies, driving 
the formation of turbulent flow at higher Reynolds num- 
bers. In parallel pipes, two simultaneous extremum prin- 
ciples can apply, depending on the flow conditions. The 
interplay of these effects is demonstrated using example 
systems. 



2. ENTROPY PRODUCTION EXTREMUM 
PRINCIPLES 

Owing to considerable confusion in the literature, it is 
important to clarify the various entropy production ex- 
tremum principles and their applicability to fluid flow in 
single and parallel pipes. Firstly, just as the purpose of 
extremum principles in thermodynamics is to predict the 
equilibrium position, as a function of the intensive and/or 



extensive variable constraints [17[ , we should expect that 
any extremum principle of interest here will enable pre- 
diction of the steady state of a flow system, as a func- 
tion of constraints on the fluxes and/or thermodynamic 
forces. Extremum methods to obtain non-steady-state 
flows are therefore not examined. Further discussion of 
this point, in connection with a classification of different 
types of systems, is given elsewhere |57H59l |. Secondly, 
the entropy production in a pipe flow (or electrical) net- 
work - although calculated using ([T]) - is of quite differ- 
ent character to that implied by (dJ, since the index j 
now refers to each conduit rather than each phenomenon 
[3 H3 i typically, only one phenomenon - such as fluid 
flow - is of interest. In consequence, extremum principles 
designed for multi-phenomenon systems need not apply 
to pipe flow networks (and vice versa). Failure to recog- 
nise this distinction has caused a great deal of confusion 
in the literature. 

It is of course recognised that the term "steady state" 
is something of a misnomer, since it refers only to a con- 
stant mean flow and not its temporal or spatial variabil- 
ity. Indeed, it has been shown that dissipative flow sys- 
tems, such as fluid turbulence and heat convection, are 
required to exhibit variability at steady state [57|, [59| . 
For consistency with other studies, "steady state" is here 
used to describe a stationary state of constant mean flow, 
to distinguish it from thermodynamic equilibrium, but 
this does not imply lack of variability of the flow. 

Examining the available extremum principles in turn: 

• Onsager's principle [23[ applies to flows in the linear 
transport regime, i.e. those in which the fluxes are lin- 
ear functions of the forces jj = LjkFk, where j and 
k are phenomenon indices and Ljk are the (constant) 
phenomenological coefficients [2l|,[22j]. The principle in- 
volves maximisation of & less the dissipation function 
<f> = ^J2ik L jk3j3ki where LJ k x is the inverse matrix 
of Ljk- The basis of the principle is obscure, but stems 
from Markovian decay of a system towards equilibrium 
(23j j . Since it requires multiple phenomena, this princi- 
ple is not considered further here. 

• Prigogine's MinEP principle l 2-5], again applicable 
only in the linear transport regime, involves minimi- 
sation of a with respect to certain fluxes j k , subject to 
a partial set of fixed forces Fk ■ The theorem is read- 
ily proved in the assumption of local equilibrium. This 
principle also requires multiple phenomena, and so is 
not considered further here. 

• The MaxEP principle of Ziegler and Paltridge is of 
quite different character, being advocated for systems 
well beyond the linear regime, subject to one or more 
fixed forces and/or fluxes. Due to the formation of 
non-linear response mechanisms, such systems can re- 
spond in many ways to the driving forces and fluxes, or 
in other words, exhibit multiple (stable) steady states. 
The MaxEP principle holds that the observed steady 
state is that which maximises the entropy production 
of the system. Aside from empirical evidence in many 
different fields (JQ), several attempts have been made 
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to derive this principle, including path-based analyses 
of transient non-equilibrium systems [6(il - [63l | , a fluctua- 
tion analysis [(34| and a relative velocity argument . 
Recently, the author [57l - l59| has given a local, con- 
ditional derivation of this principle by a probabilistic 
analysis of instantaneous fluxes through each element 
of a flow system. The analysis is identical to that used 
in thermodynamics to infer the equilibrium position, 
but adopts new entropy and free-energy-like functions 
defined for a flow system, so as to infer the steady-state 
position. Since the Ziegler-Paltridge MaxEP principle 
is considered to apply to single and multiple phenom- 
ena, it is examined here. 

A separate MinEP principle has been developed in en- 
gineering circles, especially by Bejan (§l-[To|. in which 
the entropy production of an engineered system is min- 
imised to determine the most efficient design, e.g. that 
with the highest energy efficiency. Accordingly, this is 
a design principle, which serves a different purpose to 
the predictive principles of interest here. 

Related to this discussion is an observed feature of 
many natural and engineered flow systems, described 
by Bejan as the "constructal law", in his words [Tol.l66l. 
|67I |: "For a finite-size flow system to persist in time 
... it must evolve in such a way that it provides easier 
and easier access to the currents that flow through it." 
This "law" appears to be connected in some manner to 
an entropy production extremum principle, since the 
system minimises its resistance for a given flow rate 
and/or maximises its flow rate for a given resistance. 
Although important, evolvable flow systems are not ex- 
amined further in this study. 

A separate class of MinEP principles has been devel- 
oped in finite-time thermodynamics, to calculate the 
fundamental limit of efficiency of a process in which a 
thermodynamic system is moved from one equilibrium 
position to another, along a specified path at specified 
rates (H-fll]. The geodesic is used to determine the 
most efficient path. This method has been generalised 
to steady-state flow systems 73]. Of course, it is always 
possible to do worse (produce more entropy) than these 
limits. Since this method relates to transitions between 
equilibrium or steady states, not the prediction of the 
states themselves, it is not examined further here. 

A separate MinEP (or minimum heat production) prin- 
ciple has been _advocated for the analysis of electrical 
~ 3 



circuits 13 
works [ill 
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and more recently for fluid flow net- 
12|. This involves minimising (fTJ - with j 
now indicating the pipe index - with respect to each 
flow rate, subject to the set of continuity equations 
(conservation laws) at each node. This "pipe network 
MinEP principle" is examined in detail here. 

A separate MaxEPMod principle has recently been 
given for electrical circuit analysis f27T - l30l ]. In this 
method, (TT]) is maximised - again for pipe index j - 
with respect to the flow rates, subject to the constraint 
that the heat production equals the rate of work on the 



system at steady state. Some workers [20( have con- 
fused this approach with the Ziegler-Paltridge MaxEP 
principle, when in fact it is closer in character to the 
pipe network MinEP principle. This approach is exam- 
ined herein. 

• A set of methods for pipe network analysis have been 
proposed based on maximisation of the Shannon en- 
tropy (MaxEnt), as defined in several ways [51rl56l ]. 
Although developed more for the purpose of reliability 
engineering and optimal network design (operational 
research), these methods are commented upon here. 



3. PRINCIPLES OF FLUID FLOW IN PIPES 

We now consider one-dimensional, isothermal 1 , con- 
stant composition, steady-state flow of an incompress- 
ible fluid between sections A and B in a single pipe, 
which satisfies the followin g energ y conservation (aug- 
mented Bernoulli) equation [74l . |75| : 



za^ h a A — = z B H \- a B — + H L 

2.9 pg 2g 



(2) 



where z.;, pi and Ui are respectively the elevation, pres- 
sure and mean velocity of the flow at section i, a>i is the 
kinetic energy correction factor at section i (which cor- 
rects for non-one-dimensional effects) , p is the fluid den- 
sity, g is the acceleration due to gravity, and Hl = Pl/ pg 
is the head loss between sections A and B, equal to the 
pressure loss pl per unit specific weight of the fluid. Note 
that pl does not simply give the pressure drop pb — Pa, 
but represents the energy per volume lost due to fric- 
tion between A and B. We here consider only long, 
straight pipes of uniform circular section, of pipe diame- 
ter d, mean velocity V a = Ub = U and volumetric flow 
rate Q = ird 2 U /4. From the second law of thermodynam- 
ics, in all cases Hl > in the direction of flow. However, 
to allow either choice of flow direction, it is convenient 
to consider Hl > (head loss) for Q > and Hl < 
(head gain) for Q < 0. The head loss is given by the 
Darcy-Weisbach equation: 



H r 



2gd 



u\u\ 



8fL_ 

n 2 gd 5 



Q\Q\ 



(3) 



where / is the Darcy friction factor, a dimensionless drag 
force (= 4 times the Fanning friction factor), and L is 
the pipe length. The absolute value signs in ([3]) allow for 
reverse flow. 



1 Frictional flow causes heating, and so cannot strictly be isother- 
mal. In keeping with most fluid mechanics references, the result- 
ing temperature change in incompressible flow is here considered 
negligible [see|^, p 37]. 
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In laminar flow, the friction factor and head loss can 
be derived analytically as 0, [zH] : 



fla 



H 



L,lam 



64 
128pL 



Q 



npgd 4 

in which the Reynolds number of flow is given by: 



Re 



pdU ApQ 



irpd 



(4) 
(5) 

(G) 



(allowing Re 0) and p, is the dynamic viscosity. Eq. ([5]) 
is the well-known Hagen-Poiseuille equation. For turbu- 
lent flow, no theoretical solution is available, and so the 
friction factor is usually correlated using the following 
semi-empirical relations, respectively for flow in a smooth 
pipe, in a fully rough pipe or in general [74], [75j : 



1 



= 2.01og 10 (|i?e|ViQ-0.8 



fb = -2.01o gl0 (i) + 1.14 



(7) 
(8) 



2.0 



turb 



log 10 (~) 



1.14 



d 

2.01og 10 



9.28 



\Re 



(9) 



where e is the equivalent sand roughness, a measure of the 
pipe wall roughness. Eqs. (H])-© can be solved in con- 
junction with ((3]) to give the turbulent head loss. How- 
ever, for engineering purposes, they are usually applied 
graphically using a plot of / against |i?e| (the "Moody 
diagram"). This diagram, constructed from eqs. (j4|) and 
©-(|H]), is shown in Figure [T] As shown, the flow remains 
laminar up to approximately |i?e| ~ 2100, at which a 
transition takes place to a flow of fluctuating laminar 
and turbulent character. Typically, by |i?e| « 4000 the 
flow has become entirely turbulent, although laminar flow 
can be prolonged to |i?e| > 10000 in careful, vibration- 
free experiments 0,0|- Once the flow becomes turbulent, 
for non-zero e it follows one of the curves f(\Re\,e) (|9|) 
with increasing |i?e| until it reaches a condition of "fully 
developed turbulent flow" ([5]), whereupon / becomes in- 
dependent of | i?e | and is a function only of e/d. 

Flow in a smooth pipe (e/d —> 0) follows the curve 
given by 0. This is often correlated using the simpler 
Blasius relation, valid over 4000 < \Re\ < 10 5 0,[li| : 



Bla 



0.316 

|i?e|°- 25 



which in © gives: 

Hh,Blas 



1.7&7p 025 L 

^1.75^4.75^0.25 



Q\Q\ 



0.75 



(10) 



(11) 
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FIG. 1: Plots of f(Re) for laminar flow ((4} and turbulent flow 
in smooth ([7]) and rough ([9]) pipes ("Moody diagram"). The 
Blasius correlation (|10[l is also shown. 



Eq. (TTU1) is also shown in Figure [TJ In hydraulics, the 
laminar, Blasius and fully developed turbulent head loss 
curves are often expressed using the power law: 



H L , a = XQ\Q\" 



(12) 



where the coefficient X is a function only of the pipe and 
fluid properties, whilst 1 < a < 2 expresses the transition 
between laminar (a = 1) and fully developed turbulent 
flow (a = 2). 



4. ENTROPY PRODUCTION IN A SINGLE 
PIPE 

In all the above cases, the rate of entropy production 
a per unit specific weight of fluid for incompressible, 
steady-state flow in a single pipe, at constant temper- 
ature T (equal to the power loss Pl per unit specific 
weight), is given by |8T-[l0l]: 



Pl 
99 



aT 
99 



\H L Q\ > 



(13) 



The dimensions of (|13l) suggest use of the following di- 
mensionless group (normalised rate of entropy produc- 
tion or power loss per unit length of pipe) : 



T = 



4Pl 



4rjT 



TidpgL ixd[igL 



Hr 



-Re 



1 |i?el 3 



where the Galileo number Ga — p 2 gd 3 /p 2 , which ex- 
presses the ratio of gravity and inertial to viscous forces, 
depends only on the pipe and fluid properties [7^, l8lj ]. 
To the author's knowledge, the composite group T (fbij) 
has not been suggested previously. Several workers 
[llT [l2l [75j use the group fRe 2 or its square root; how- 
ever, this only correlates with the dimensionless head loss 
Hl/L — i/ Re\Re\/Ga, not the entropy production. 
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FIG. 2: Plots of (a) TGa and (b) T/T iam against Re for the 
same curves as in Figure [T] calculated using (|14|l . The insets 
show the transition region. 



Plots of the normalised entropy production TGa 
against |i?e| for the various curves for laminar and tur- 
bulent flow are shown in Figure Ufa). Plots of the ratio 
T/T; am against |i?e| are given in FigureHJb). As evident, 
for |i?e| < 1034 the entropy production by laminar flow 
in a smooth pipe exceeds that of turbulent flow, but the 
latter increases more rapidly with \Re\; for |i?e| = 10 4 , 
10 6 and 10 8 , the respective ratios are T/T/ am = 4.82, 
182.0 and 9282. The known onset of the transition re- 
gion at |i?e| w 2100 is of the order of (a little higher 
than) the crossover point (see later discussion). Turbu- 
lent flow in rough pipes follows similar trends, with even 
higher entropy production than a smooth pipe. Yet, for 
all Reynolds numbers, the laminar solution (lU)-© IS a 
valid steady-state solution, being the analytical solution 
of the Navier-Stokes equations. The laminar and tur- 
bulent curves in Figures dta)-(b) therefore represent two 
alternative steady-state solutions, from which the system 
appears - for the most part - to select the one with the 
highest entropy production. Thus although not consti- 
tuting a proof, Figures [UJa)-(b) provide strong evidence 
for the action of the Ziegler-Paltridge MaxEP principle 



in providing a driving force for the selection of lam- 
inar or turbulent flow. 

The above argument was given by Paulus Jr [ll], [l2| 
and Martyushev [l3| , but presented graphically in terms 
of the friction factor instead of the entropy production. 
Note that the trends of the entropy production curves 
in Figure [DJa) are remarkably similar to those of other 
"complex systems" which undergo a critical transition 
between linear and non-linear dissipation mechanisms. 
These include (i) heat transfer through a convectivc fluid, 
for which a dramatically increases with Rayleigh num- 
ber at the transition from heat diffusion to convection 
[80l ] : and (ii) chemical degradation in an ecosystem, for 
which & dramatically increases with chemical affinity at 
the transition from chemical to biological metabolic pro- 
cesses [13] ■ Clearly, the Ziegler-Paltridge MaxEP concept 
unites many quite different disciplines (see [57|). 

Examining the discrepancy between the typical onset 
of turbulent flow at |i?e| « 2100 and the crossover point 
of |i?e| = 1034 in Figures (Ha)-(b) (with lower values in 
rough pipes), two explanations are forthcoming. Firstly, 
the discrepancy may result from the use of empirical re- 
lations for the turbulent friction factor 0-©, which are 
not fitted to nor intended to be applied to the transition 
region. Secondly, the discrepancy may be due to the for- 
mation of "metastable" laminar flow at higher Reynolds 
numbers than would otherwise be expected from the en- 
tropy production. This phenomenon is well known, as 
demonstrated by the sudden transition to turbulence dur- 
ing flow in a pipe, in contrast to the more gradual tran- 
sition seen in flow in packed beds and many external and 
boundary layer flows @, 0, 0, [lH, [8l|. Also, as men- 
tioned earlier, laminar flow can persist at higher |i?e| in 
careful experiments, even above 10000, but is unstable to 
vibrational disturbance [6, 7] . This situation is somewhat 
analogous to the existence of a metastable supersaturated 
solution of a salt, despite a free energy gradient which 
should drive its precipitation; in such cases, perturbation 
of the system (such as tapping the container wall) can suf- 
fice to drive it to equilibrium. No-one discards thermo- 
dynamics because of this exceptional phenomenon. The 
relevant question here is not whether metastable laminar 
flow can occur at high |i?e|, but what is the minimum 
critical |i?e| for the onset of turbulence? Some research 
in this direction has been undertaken by Benhamou et al. 
[82| [see also , who reported significantly lower criti- 
cal |i?e| of around 1200-1500 in a horizontally oscillating 
pipe, and by Li et al. [83j, who found the critical |i?e| 
can drop to 1600-2000 due to pipe entrance effects. 

For completeness, we also consider the (somewhat ab- 
surd) variation of the pipe entropy production with re- 
spect to the flow rate Q. Taking ([TB")) as the objective 
function F, using the power-law relation (fT2")l for the head 
loss, differentiation gives: 

pi -p 

— = (a + l)XQ\Qr 1 = (a + l)H L (15) 
This is non-zero unless Q = 0, hence the steady state is 
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Q 2 J 2 

FIG. 3: A simple parallel pipe network with fixed total flow 
rate A. 

not an extremum in the entropy production with respect 
to Q. This does not affect the above-mentioned MaxEP 
principle, which applies with respect to the selection of 
flow regime and is unrelated to the flow rate. However, 
(fT5| indicates that the pipe network MinEP principle - 
which involves the variation with respect to flow rates - 
does not apply to fluid flow in a single pipe. The validity 
of the MaxEPMod principle of Zupanovic and co-workers 
for flow in a single pipe is deferred to ^5.21 

5. ENTROPY PRODUCTION IN PARALLEL 
PIPES 

5.1. Fixed Flow Rate Systems and the Pipe 
Network MinEP Principle 

We now turn attention to the simple pipe geometry of 
Figure |31 in which the flow Q is forced to subdivide be- 
tween two parallel pipes, of flow rates Q\ and Q 2 , which 
then rejoin. The pipes have different flow characteristics. 
We here neglect head losses associated with the pipe junc- 
tions and fittings. Given the pipe and fluid properties, 
how will the flow divide between Qi and Q 2 1 In fact, a 
closed-form solution exists. Firstly, the flows are subject 
to continuity (conservation of fluid) equations at each 
node, in this case identical at both nodes [74j : 

Q-Qi-Q 2 = (16) 

This is equivalent to Kirchhoff 's current law in electrical 
circuit analysis, taking as analogues the current Ij and 
flow rate Qj [ll, . Secondly, the system is considered 
to satisfy the condition that at steady state, there should 
be no circular flow around any loop, or alternatively, that 
the head losses along all paths between two nodes must 
be the same [74j : 

H L i -H L2 ={) (17) 

This is a non-linear variant of KirchhofFs voltage law in 
electrical circuit analysis [2^ [ill ; taking the voltage drop 
Vj and head loss Hlj as analogues, consistent with the 
notion that each node should be at a constant steady- 
state piezometric head (potential). Substituting ([3]) for 
the head losses in ([T7| . using the friction factor (J4j) , ([Tj) , 
J5]) or © for the applicable flow regime, yields two equa- 
tions in the two unknowns Q\ and Q 2 . For example, 



using the power law relation (fl2|) . (|17l) reduces to 

XiQ^Q^- 1 -X 2 Q 2 \Q 2 \ a - 1 =0 (18) 

Knowing a and X, for j 6 {1,2}, (fT6|) and (fT8]) can be 
solved numerically for Q\ and Q 2 - If the flow regime 
is not a priori known, the solution requires an iterative 
selection scheme based on calculation of Rej and Xj. 

Although a closed-form solution exists, its solution 
becomes inconvenient for larger pipe networks, requir- 
ing continuity equations for all nodes and head balance 
equations for all loops. Usually, the problem is over- 
constrained (more equations than unknowns), causing 
difficulties with numerical solution. However, there exists 
a quite different approach, which appears to have been 
almost completely overlooked in the enginee ring litera- 
ture, except for two studies by Paulus Jr [ill, EH- A lin- 
ear variant of the method, for electrical circuit analysis, 
is slightly better known 0, 131rf33j . The method involves 
constrained minimisation of the total entropy production 
at constant T (equal to the total power loss) within the 
pipe system: 

99 99 j~[ 

with respect to the flow rates, subject to the node conti- 
nuity constraints: 

J 

Qn + J2^n = 0, n = l,...,N (20) 
j=l 

where j is the pipe index, J is the number of pipes, n 
is the node index, N is the number of nodes, Qj n is the 
flow rate from pipe j into node n and Q n is the external 
inflow into node n. 

It is worth examining the validity of this method for 
parallel pipe networks. For the simple example in Figure 
131 assuming power law behaviour (|12p with constant a, 
proof is relatively straightforward. From flTT)l) : 

^ = — = \H Ll Q 1 \ + \H L2 Q 2 \ (21) 
99 99 

Substituting for the head losses (|T2"j) and Q 2 (fTl)]) gives 
the constrained objective function: 

F = X 1 |Q 1 | a+1 + X 2 \Q - Qx| a+1 (22) 

Differentiation yields: 

BF 

— = (a + l^QrlQrr 1 - X 2 (Q - Q X )\Q - Qx^" 1 ] 
0Q1 

= (a + l^QrlQrr 1 - X^Q^ 1 ] (23) 
= (a+l)[H L1 -H L2 ] 

where (fT2"]) and (fTl))) are again used. For steady state flow, 
the last form of (|2"5)) must be zero, hence the steady state 
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is an extremum. Differentiating a second time gives: 



d 2 F 



aia + l^lQ^- 1 



x 2 \Q-Qi\ a - 1 ] 



aia + l^lQ^- 1 +X 2 \Q 2 \ a - 1 } 



a(a + 1) 



a— 1 

Hli , a li 



again using (fT2")l and (ITC)) . Since the terms Hi,j/Qj are 
always positive at steady state irrespective of flow direc- 
tion, d 2 F/dQ\ > for a > 1, and so the extremum is 
indeed a minimum. A related proof, only for the lami- 
nar case, was given by Paulus Jr [llj. Since this method 
uses the node continuity constraint (fll)]) . and from (|2"3"1) . 
recovers the loop head constraint (flTl) . it is identical to 
the closed-form solution. 

The above proof can be extended to any number J > 2 
of parallel pipes with common a between two nodes, in 
which case one obtains J — 1 derivatives with respect to 
Qk, each of which is a local minimum. The steady state 
position is therefore a global minimum. The MinEP and 
closed- form solutions are again equivalent. 

Difficulties arise, however, with the pipe network 
MinEP principle for other forms of the friction factor 
f(Re). Consider the case of two parallel pipes (Figure[3]) 
with constant but different power-law coefficients a\ and 
a 2 in (HI]). From (16]) and (fl9|): 



It is important to clarify the differences between the 
)ipe network MinEP principle and that of Prigogine 0- 
5]. Firstly, as noted earlier, the former involves the con- 
strained minimisation of a sum of entropy production 
(24) terms for J > 2 pipes, within which the flow rates Qj are 
not independent, whereas Prigogine's theorem involves a 
sum of entropy production terms for J > 2 independent 
physical phenomena. Secondly, the present theorem ap- 
plies to both linear (a = 1) and non-linear (1 < a < 2) 
functional relationships H^jlQj), albeit only for the re- 
stricted case of systems with common a. This is in con- 
trast with Prigogine's method, which requires linear (On- 
sager) relations between the forces and fluxes. However, 
the present theorem is in some sense allied to Prigogine's, 
in that both emer ge - in some circumstances - as a sim- 
ple (almost trivial [33j |) consequence of stationarity of the 
steady state compared to non-steady-state choices of the 
flow rates. As such, the pipe network MinEP principle - 
when it applies - relates only to the dynamic stability of 
the steady state, and does not have any bearing on the 
selection of the steady state if there exist multiple steady 
state solutions. 



F = Xi|Qi| Ql+1 + X 2 \Q - Q x 



Q2 + 1 



from which 
dF 



dQi 



= (ai + l)X 1 Q 1 \Q 1 



ai-l 



(25) 



(26) 



- (02 + 1)A 2 (Q - Q X )\Q - Qxp- 1 
= ( ai + l)X 1 Q 1 \Q 1 \ a ^ 1 - (a 2 + l)X 2 Q 2 \Q 2 \ a "- 1 
= ( 0l + l)H Ll - (a 2 + l)H L2 

This is certainly not zero unless (trivially) a\ — a 2 or the 
pipes carry flows with Hl\ ^ H^ 2l an unstable situation. 
In consequence, the pipe network MinEP principle does 
not apply to parallel pipes of different aj. Indeed, as 
proven in Appendix [A] it also fails for parallel pipe net- 
works with non-power-law /(-Re) relations such as those 
in ([7]) or ([9]). Thus although it is quite appealing to con- 
sider the pipe network MinEP principle as a fundamen- 
tal principle of fluid mechanics, as advocated by Paulus 
Jr [ll], [ljj, it is not. A similar conclusion was reached 
by Landauer 32j with regards to electrical networks, for 
which the MinEP principle applies to circuits contain- 
ing resistors, but not those with inductances. A further 
difficulty arises for variable temperature flows, which do 
not satisfy any MinEP principle, a criticism also made by 
several authors |lj§ [32], |33j in reference to electrical net- 
works. The term "minimum entropy production princi- 
ple" is therefore rather misleading, since the principle un- 
covered here is more accurately described as a "minimum 
power loss" principle, applicable only to the restricted set 
of isothermal flows in constant a pipe networks. 



5.2. Fixed Head Systems and the MaxEPMod 
Principle 

We now consider a new method for electrical cir- 
cuit analysis given by Zupanovic and co-workers (27l - [30| , 
which is here extended to a parallel pipe network. A 
slightly different flow configuration is required, as repre- 
sented in Figure HJa), in which the system is constrained 
by a fixed change in head AH (assuming infinite reser- 
voirs) rather than a fixed total flow rate Q. The analo- 
gous electrical network is shown in Figure|4][b). From the 
head loop principle, the system yields two loop equations: 



AH = H 



J = 1,2 



(27) 



which can be solved for Q\ and Q 2 . The MaxEPMod 
method involves maximisation of the entropy production 
multiplied by constant T (equivalent to the heat pro- 
duction q) (f2"Tj) , subject to the constraint that the heat 
production is equal to the work expenditure w: 



w = pg\AHQ\= pg\AH\\Q 1 + Q 2 



(28) 



Applying the calculus of variations, assuming power-law 
pipes (fTSj) of constant a, gives the Lagrangian: 



G = JL + A( - _A H(1 _ A) A + A - 

99 99 99 99 99 



(29) 



= (1 - A)(|ff L iQi| + \H L2 Q 2 \) + X\AH\ |Qi + Q 2 \ 
= (1 - A)(X 1 |Q 1 | a+1 + X 2 \Q 2 \ a+1 ) + X\AH\ \Q X + Q 2 \ 

where A is the Lagrangian multiplier. Setting the varia- 
tion with respect to each Qj to zero gives: 
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(b) 



Q 2 ^ 2 



I — YA — I 

/ 2 ,f? 2 



FIG. 4: Network diagrams for the MaxEPMod method: (a) 
parallel pipes between reservoirs, and (b) analogous electrical 
circuit. 



dG 



{1- \)(a + l)X ] Q ] \Q ] \ a ~ 1 



A|Ai7|sign(Q 1 + Q 2 ) =0 (30) 



which, from (fT2|) . gives: 



(a + l)fT 



(a + l)H Lj - |Aff|sign(Q)' 



Y? (31) 



For Q ^ 0, the externum always coincides with the loop 
equations AH = \AH\ sign(Q) = H^j^j, if we set A 
a + 1)/ a. This reduces to the value of A = 2 given 
previously for electrical networks (a = 1) [27], HH. Since 
the method recovers the two loop constraints - the node 
continuity now being irrelevant - it is identical to the 
closed-form solution. Taking the second derivative gives, 
for Q j £ 0: 



d 2 G 



dQ 



2 = (l-X)a(a + l)X J \Q J \ a ' 1 



-(o + l)A>|Q 3 -| e 



(a + l)H Lj 



(32) 



Again H Lj /Qj > 0,Vj and so d 2 G/dQ\ < for a > 
1, hence the extremum is a maximum. By symmetry, 
the above analysis can be generalised to any number of 
parallel pipes, and indeed, even to a single pipe. 

Difficulties again emerge, however, in other situations. 
For pipes with different a.j, one again obtains (|3"0"j) - (f3"Tj) 
but expressed in terms of Oj instead of a. No com- 
mon value of A exists which can satisfy all loop con- 
straints. More broadly, as proven in Appendix [B] it 
also fails for parallel pipe networks with non-power-law 
f(Re) relations. Furthermore, just like the pipe network 



MinEP principle, the MaxEPMod approach fails for vari- 
able temperature systems, and so is more accurately de- 
scribed as a heat-work extremum principle 2 . The MinEP 
and MaxEPMod principles therefore have many features 
in common, and so can be regarded as kindred principles 
differing primarily in their choice of flow conditions. 



5.3. MaxEnt Methods 

As mentioned earlier, several methods for pipe network 
analysis have been proposed based on maximisation of 
the Shannon entropy (MaxEnt) [5ll - l56l |. Of these, most 
have been developed for reliability or sensitivity analy- 
sis, e.g. the method of Awumah and co-authors [5l|, in 
which the pipe network of highest entropy is considered 
the most reliable. Tanyimboh and Templeman [52[ pro- 
pose three separate MaxEnt methods for the analysis of 
separate types of pipe networks, with various definitions 
of probabilities and entropies. Unfortunately, all three 
methods permit only one connection between each pair 
of nodes, and therefore cannot be applied to the parallel 
pipe problems considered here. Furthermore, none of the 
methods incorporate frictional resistances for each pipe 
(such as fj or Xj) and so their results assume equivalence 
of all pipes, a very unrealistic situation. 



5.4. MaxEP Principle 

For parallel pipe flow, we must not forget the action of 
the quite distinct Ziegler-Paltridge MaxEP principle (<J2|). 
The interplay between this extremum principle and the 
MinEP or MaxEPMod principles is examined in the next 
section, by the analysis of examples for constant flow rate 
or constant head conditions. 



6. EXAMPLES 

We now illustrate the various roles of the Ziegler- 
Paltridge MaxEP, pipe network MinEP and Zupanovic 
MaxEPMod principles in parallel pipe networks by two 
examples. Example 1 is of two parallel pipes with fixed 
total flow rate Q (Figure [3]), in which the pipes can expe- 
rience either laminar flow (U])-© or turbulent flow rep- 
resented by the Blasius relation (pTTj) - (fTT j) . Since these 
relations conform to power-law behaviour (|12[) . the pipe 
network MinEP principle applies; we also examine the 
Ziegler-Paltridge MaxEP principle. The system is as- 
sumed to carry water, with properties p — 1000 kg m~ 3 , 



2 Note from the first line of ((29), knowing A > 1, the MaxEPMod 
principle could equally have been posed as a MinEP principle, 
subject to a constant work constraint, or a maximum work prin- 
ciple, subject to a constant heat production constraint. 
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/i = 10~ 3 Pa s, g = 9.81 m s -2 and with both pipes of the 
same diameter di = d 2 — 0.05 m; the difference between 
pipe resistances is due to their different lengths, L\ — 50 
m and L 2 — 150 m. In consequence, the Reynolds num- 
ber of each pipe Rej is proportional only to its flow rate 
Qj and independent of diameter dj (j6]), whence Q can 
be expressed in dimensionless form in terms of the com- 
posite Reynolds number Re to t = ^pQ/n^id = Re\ + Re 2 . 
In each case, the laminar and turbulent steady-state so- 
lutions over a range of Retot were separately computed 
by minimisation of the normalised entropy production 
(|21|). subject to the continuity constraint (TTB)) . Compu- 
tations were conducted at 15-digit precision using the 
Minimize command in Maple 12 [? ], which invokes a 
linear or nonlinear optimisation routine as needed; in all 
cases, convergence was attained within 4 iterations for 
both laminar and turbulent flow. Cross-checking of each 
solution against its closed-form counterpart (jT^j) - (p~7)) in- 
dicated flow rate errors of < 10~ 19 m 3 s _1 for laminar 
flow and < 10~ 9 m 3 s~ x for the Blasius relation. 

The calculated laminar and turbulent steady-state so- 
lutions for Example 1 are illustrated as solid lines in 
the three-dimensional plot of normalised entropy pro- 
duction against Retot and Re\ in Figure Eta). As ex- 
pected, the laminar curve is of higher entropy produc- 
tion at lower Reynolds numbers, but is overtaken by the 
turbulent curve at higher Reynolds numbers 3 , consistent 
with the Ziegler-Paltridge MaxEP principle. Superim- 
posed on this plot are several subsidiary sets of curves 
for specific values of Re to t, which illustrate the locus of 
non-steady-state solutions for laminar or turbulent flow, 
as obtained by simultaneous solution of (f2~Tj) and the con- 
tinuity constraint (fT6|) for given values of Retot and Re\. 
As evident, in every case - both for laminar and turbu- 
lent flow - the steady-state solution lies at the minimum 
of this locus. Figure [Sfa) therefore illustrates the simul- 
taneous action, in parallel pipe systems of fixed Retot 
(or Q), of (i) the pipe network MinEP principle, which 
discriminates the steady-state position from the set of 
non-steady-state solutions for each flow regime and spec- 
ified Retot , and (ii) the Ziegler-Paltridge MaxEP princi- 
ple, which provides a means to determine the flow regime 
at a specified Retot- 

The next Example 2 also involves two parallel pipes, 
but with a fixed total head AH (Figure Sta)) rather than 
fixed Q. The fluid and pipe properties are otherwise iden- 
tical to those in Example 1. We again consider laminar 
©-© or Blasius (|TD ]) -(|TT 1) flow, both of power law form, 



3 Since the lines do not intersect, for equality of a and satisfaction 
of the continuity and head loss constraints l|16ll - l|17|l . the system 
will exhibit a transition region over the range Retot = 1585 to 
3417, with turbulent flow in pipe 1 and laminar flow in pipe 2. 
In each case the transition to turbulence is predicted to occur 
in each pipe at Rej = 1189. If the smooth curve J7} is used 
instead of the Blasius curve, the predicted transition values shift 




locus of 

non-steady-state 
turbulent points 
locus of 

non-steady-state 
laminar points 



locus of 

non-steady-state 
laminar points 



locus of 

non-steady-state 
turbulent points 



AH (10-3 m) 



FIG. 5: Plots of normalised entropy production in parallel 
pipes for (a) constant | .Retot | (Example 1), showing action of 
Ziegler-Paltridge MaxEP and pipe network MinEP principles, 
and (b) constant AH (Example 2), showing action of inverted 
Ziegler-Paltridge MaxEP and Zupanovic MaxEPMod princi- 
ples. For details of calculations, see >J6] 



hence the MaxEPMod principle applies. In this case, the 
laminar and turbulent steady-state solutions over a range 
of AH were computed by maximisation of the normalised 
entropy production (or heat generation) (|2ip. subject to 
the heat-work constraint, which from (|2f I) and ([28]) can 
be expressed in the form: 



*i (<0 |Qi| o+1 



-X { 2 a) \Q 2 \ a+1 = \AH\\Q 1 + Q 2 \ (33) 



to Retot = 1379 to 3039 and Re, 



1034. 



where a = 1 and 1.75 respectively for the laminar and 
Blasius relations, and Xj are constants for each spec- 
ified pipe and flow regime. The Maximize routine in 
Maple f2 with 15 digit precision was used for compu- 
tations; in this system, convergence required 30 to 50 
iterations, with an error of < 10~ 7 m 3 s _1 . The cal- 
culated laminar and turbulent steady-state solutions are 
illustrated as solid lines in Figure [S^b), now as a function 
of AH and Re\. In contrast to Example 1, the lami- 
nar curve is of lower entropy production than the Blasius 
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curve at lower AH, but overtakes it 4 at higher AH. The 
analysis therefore suggests the action of a MinEP princi- 
ple for the selection of flow regime, quite opposite to the 
Ziegler-Paltridge MaxEP principle. However, under con- 
stant AH conditions, the flow rates are not comparable; 
for example, AH = 0.003 m produces a total laminar 
flow of Q = 1.20 x 10~ 4 m 3 s" 1 (Re to t = 3066) or a tur- 
bulent flow of only 1.04 x 10" 4 m 3 s" 1 (Re tot = 2659). 
The analysis therefore confirms the argument of Paulus 
Jr and co-authors [Tl|, [l2| that under constant head (or 
pressure) conditions, the flow regime can be selected by a 
MinEP principle. Since the transition values are identi- 
cal to those given for Example 1 by the Ziegler-Paltridge 
MaxEP principle, it is a simple inversion of the latter, 
not an entirely new principle. 

As in the previous example, non-steady-state curves for 
specific values of AH are also superimposed on the plot, 
calculated by simultaneous solution of ([2"TT) and the heat- 
work constraint (|33|) for given values of AH and Q\. As 
evident, the steady-state solution lies at the maximum of 
this locus, both in laminar and turbulent flow. A curious 
feature of this system is that the subsidiary curves are 
one-sided, i.e., the non-steady-state flow rate Qi ( or Rei) 
cannot exceed its value at steady state. Examining (|33)) . 
we see that this describes a piecewise convex function 
for AH(Qi,Q2) > (convex within each quadrant in 
Qi and Q 2 , centred on the origin), which for constant 
AH does not permit any real solution for Qi higher than 
the steady-state value. Figure [SJb) thus illustrates the 
simultaneous action, in parallel pipe systems of constant 
AH, of (i) the MaxEPMod principle of Zupanovic and 
co-workers, which discriminates the steady-state position 
from the set of non-steady-state solutions for each flow 
regime and specified AH, and (ii) a strict inversion of 
the Ziegler-Paltridge MaxEP principle, which provides a 
method to determine the flow regime at a specified AH. 



"complex systems" such as heat convection [8(| and bio- 
logical metabolism [47j . suggesting wide applicability of 
this principle. Secondly, flow in two or more dissimilar 
parallel pipes subject to a constant total flow rate can 
simultaneously satisfy two principles: (i) a pipe network 
MinEP principle with respect to the choice of flow rates, 
and (ii) the Ziegler-Paltridge MaxEP principle with re- 
spect to the choice of flow regime. The former princi- 
ple - different to but allied to that of Prigogine - arises 
as a simple consequence of the stability of the steady 
state compared to non-steady-state flows; it is proven 
for isothermal laminar and turbulent flows in parallel 
pipes with a constant power law exponent, but is oth- 
erwise proven to be invalid. The latter principle appears 
to be more fundamental. Flow in parallel pipes subject 
to a constant head condition can satisfy: (i) a modified 
maximum entropy production (MaxEPMod) principle of 
Zupanovic and co-workers with respect to the choice of 
flow rates, and (ii) an inversion of the Ziegler-Paltridge 
MaxEP principle with respect to the flow regime. The 
former is again valid only for flows with a constant power 
law exponent. 

In either type of parallel pipe system, the pipe network 
MinEP or MaxEPMod principles provide the same infor- 
mation as the known closed-form solution and, as shown, 
are not of universal validity. Such methods may therefore 
be useful for numerical solution purposes, but cannot be 
considered as fundamental principles. (Similar consider- 
ations apply to the electrical flow analogues of these prin- 
ciples.) In contrast, the Ziegler-Paltridge MaxEP prin- 
ciple or its inversion provides a criterion with which to 
choose between two known steady-state solutions; it is 
therefore of much greater importance for the analysis of 
fluid flow systems. 

Further work is required on the validity of entropy pro- 
duction extremum principles in more general pipe net- 
works, including those with series and parallel flows and 
with inflows or outflows at each node. 



7. CONCLUSIONS 

The present study examines the role of various en- 
tropy production extremum principles for the analysis of 
steady-state flow of an isothermal, incompressible fluid 
subject to friction in single or parallel pipes. Firstly, 
flow in a single pipe is shown to satisfy the Ziegler- 
Paltridge MaxEP principle, which provides a criterion 
for the choice of laminar or turbulent flow. The result- 
ing plots of entropy production against Reynolds number 
(indicating distance from zero-flow equilibrium) are re- 
markably similar to those reported for other dissipativc 



4 Again the lines do not cross, forcing formation of a transition re- 
gion over AH = 0.00155 to 0.00466 m using the Blasius relation, 
or AH = 0.00135 to 0.00405 m using the smooth relation. The 
total and pipe Reynolds numbers at the transition endpoints are 
identical to those for Example 1. 
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Appendix A: The Pipe Network MinEP Principle 
for Arbitrary /(-Re) 

Consider the case of two parallel pipes subject to a con- 
stant total flow rate Q (Figure with variable friction 
factors fj(Rej). From © and (fll))) . with the flow rates 
Qj expressed as functions of Reynolds numbers Rej , the 
entropy production per unit specific weight is: 

cVT _ V ( hL^Re^ f 2 L 2 \Re 2 \ 3 \ 

99 8p 3 . 9 V d\ + 4 ) [Ai) 
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From node continuity (|16p . Re 2 = kpQ l'K\id 2 — Re\d\l d 2 , 
giving the objective function: 



F 



nfi 3 fhL^Rdl 3 , f 2 L 2 



'<J 



4 



4 



ApQ 



TTfJ, 



Reidi 



(A2) 



Setting the derivative with respect to Re\ to zero, and 
back-substituting in terms of Re 2 and the head losses, 
gives: 



dF 
dRe[ 



= = 



7T//' 



3 r 



8p3g 



3Li/iJ2ci|J2ei| M x L 2 f 2 Re 2 \Re 2 \ 



4 



4 



d% 



dRe 2 



£i|i? ei | 3 8h d^Re^ 3 9/ 2 
d? dRei 

37T/idi 

-^—(H L1 -H L2 ) 
TT/i 3 rLi|i?e!| 3 a/i 



'3 



di£ 2 |i?e 2 | 3 9/ 2 



4 



9i?e 9 



(A3) 



The last form can only satisfy Ul\ = -Hl2 if it meets the 
condition: 



Li|flei| 3 9/i L 2 |i?e 2 | 3 9/2 



(A4) 



df 9i?ei d 3 , dRe 2 

or, substituting for Lj from ([3]) and the definition of Rey. 
H L1 R Cl dh H L2 Re 2 9/2 



/1 9i?ei / 2 dRe 2 
For iJ^i = -Hl 2 , this gives the condition: 
R ej dfj 



A = constant 



(A5) 



(A6) 



Eq. 



fj dRej 

6| is a differential equation with general solution: 



fi = CRej 



(A7) 



where C=constant. To preserve fj € M + for Rej < 
and C > we can, without loss of utility, restrict 
to the form: 

fj = C\Rej\ A 



(A8) 



which also satisfies (|A6|) . This proves that in a paral- 
lel pipe network, the MinEP principle conforms to the 
closed-form solution only for pipes with power-law be- 
haviour (|12p with a common exponent A — a — 2. A 
similar result ensues for any number J > 2 of parallel 
pipes. 



Appendix B: The MaxEPMod Principle for 
Arbitrary f(Re) 



Consider two parallel pipes subject to a constant head 
(Figure B|a)), with variable friction factors fj(Rej). 
From (|2T))) and the definition of Rej, the objective func- 
tion is: 



G 



+ 



(1 - \)irp 3 ( frL^Re^ 

8 P 3 g 

\\AH\irp 



f2L 2 \Re 2 \ 3 



4 



4 



(Bl) 



Ap 



R ei \dx + \Re 2 \d 2 ) 



Setting dG/dRej = and solving for A gives: 



p?LjRe*j(zf j + Rej^-) 



p?LjRe)(zfj + Rej A) _ 2d 3 p^g\AH 



L l{ :^f j + Re 3 ^ 



L] (3Re]f ]+ Re 3 ^-)-2\AH\G 



Rej dfj 
fj dRej 



Re df 3 \AH\ sign(fej) 



(B2) 



(B3) 



(B4) 



fj dRej 



H 



Lj 



where all dependent extraneous variables are eliminated 
by successively substituting for the system properties us- 
ing the Galileo number (Q, and for fj and then Lj in 
terms of the head loss, H^j = ^fjLjRej\Rej\/Ga,j (see 
(O). For the MaxEPMod principle to be valid, firstly, 
it must be true that \AH\ sign(Rej) = HLj,Vj, and sec- 
ondly, it must be possible to choose a common A for all 
pipes. Equating A for two parallel pipes then gives the 
condition: 



Re j d fj 
fj dRej 



A = constant 



(B5) 



This is the same differential equation as (|A6|) , again yield- 
ing a power-law solution (|A7|) - (|A8|) . The same result is 
obtained for any number of parallel pipes. This proves 
that in a parallel pipe network, the MaxEPMod princi- 
ple is valid only for flow in power-law pipes of common 
exponent a. 
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